The differential sum rule for the relaxation rate in the cuprates 
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Motivated by recent experiments by Basov et al, we study the differential sum rule for the effective 
scattering rate l/r(w). We show that in a dirty BCS superconductor, the area under 1/t(cj) does 
not change between the normal and the superconducting states. For magnetically mediated pairing, 
similar result holds between T < T c and T > T c , while in the pseudogap phase, 1/t(uj) is just 
suppressed compared to 1/t(uj) in the normal state. We argue that this violation of the differential 
sum rule is due to the fact that low-energy fermions in the pseudogap phase are still incoherent. 



The analysis of the optical sum rules in condensed 
matter systems is a valuable tool that helps one to un- 
derstand the key physics and relevant energy scales in 
the problem |Ej. The focus of this report are recent 
experimental results M for the effective relaxation rate 
&[1/<j( 
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i<j2{ui) is the optical conductivity, 

plasma frequency, and n is the density of particles. The 
data analysis for optimally doped Y Ba.2Cu3Os.Qb [[§ and 
TliB a^C uO§+ x [Q and underdoped Y Ba2CunO% and 
Bi2Sr2CaCu20s+ x M revealed an approximate differ- 
ential sum rule for t (uj) between T >T C and T < T c : 
although J divT^ 1 (uj) does not converge, it changes very 
little when the system enters into the superconducting 
state. This differential sum rule, however, is not satisfied 
between the normal and the pseudogap phases - 1/t(uj) 
in the pseudogap phase appears to be just suppressed. 

Quite generally, exact sum rules are related to con- 
servation laws. The /—sum rule for the optical con- 
ductivity states that at a given density of particles, the 
total absorbing power of the solid characterized by o\ 
does not depend on the details of the interactions and 
is determined only by the total number of particles in 
the system J5J . The total absorption power is given by 
Io° (J i( UJ ) Applying the Kubo formula that relates cr(u>) 
with the full retarded current-current correlator H(uj): 
(t(ui) = (cj 2 ;/47t) n(w)/(— iui + 6), separating the fre- 
quency integral into the integral over infinitesimally small 
u> and the rest and using the Kramers-Kronig relation for 
II(w) — 1 that vanishes at the highest frequencies, we ob- 
tain J °° <7i(w) = w 2 ;/8, independent of the form of n(w). 

Is there an analogous sum rule for 1/t(lo)1 Using 
1/t(u>) = — Im[Lu 2 /T1(uj)}/uj and applying the Kramers- 
Kronig relation, we find 
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The constant C again has to be chosen such that 
uj 2 /H(ui) + C vanishes at u> — > 00. However, C turns 
out to be infinite as at high frequencies II (u> ) rs 1, and 
lo 2 /Xl{u) diverges. This divergence implies that there is 
no conservation law associated with the relaxation rate 
and hence no sum rule for 1/t{u>). 

Still, one can try to see whether one can get useful 
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FIG. 1: The results for (Ti(u) and l/r(w) (in arbitrary units) 
for a BCS superconductor, to first order in 1/rA (clean limit). 
In the dirty limit, the jump in the conductivity at ui — 2A 
is much smaller. The inserts show Iit(uj) = 8/Wpj fg dxo~\(x) 
axi&I T {oj) — (t/2A) f Q dx(l/r(x)— 1/r). The arrow indicates 
the value of / CT (oo) and . The differential sum rule for l/r(u) 
is satisfied if I T (oo) — 0. 



information by comparing 1/t(lu) for two different sys- 
tem parameters, e.g. temperatures, which do not af- 
fect the system behavior at high frequencies. Suppose 
momentarily that 1/t(w,Ti) — 1/t{uj,T2) vanishes at 
high frequencies. The Kramers-Kronig relation is then 
fully satisfied and using (|l|) we immediately find that 
J du(T— 1 (w,Ti) - t- 1 (uj,T 2 )) = 0, i.e., the area un- 
der 1/t(u>) does not change with T. The conservation 
of the total area under 1/t(lu) would create a valuable 
tool to study the evolution of the spectral weight be- 
tween, e.g., the normal and the superconducting states. 
This new differential sum rule, however, is not associ- 
ated with a conservation law and therefore is not guar- 
anteed to be satisfied - only explicit calculations can 
determine whether or not the temperature dependence 
in n(w, T) is weak enough to ensure the convergence 
of 1/t(w,7i) — 1/t(o;, T2). In this communication we 
study under which conditions the differential sum rule 
for 1/t(u>) actually holds, and at which frequencies it is 
exhausted. We consider the magnetic scenario for the 
pairing in the cuprates, and argue that although the dif- 
ferential sum rule is approximately satisfied for all T, it 
is exhausted at energies comparable to the pairing gap 
only if there is a strong feedback effect from the pairing 
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on the fermionic propagator - otherwise, the sum rule is 
exhausted at much larger frequencies. We associate the 
hrst regime with T < T c , and the second one with the 
pseudogap phase. 

To put our analysis of the spin mediated pairing into 
perspective we first analyze the situation in a dirty BSC 
s— wave superconductor at T = 0, when the pairing 
causes a strong feedback on the fermionic propagator, 
and in the toy model where there is no feedback from the 
pairing on the fermionic self-energy. 

The theory of a dirty superconductor is well devel- 
oped H 0. In the normal state, the inelastic scatter- 
ing by impurities yields a retarded fermionic self-energy 
= i/2r, In a superconducting state, this self-energy 
is modified due to a feedback from superconductivity and 
takes the form S(w) = (i/2r)u)/\/u) 2 — A 2 where A is the 
superconducting gap Q. Substituting these forms into 
the current-current polarization bubble and performing 
the momentum integration we obtain 
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where ft± = ft ± lo/2. In the normal state, this reduces 
to a conventional Drude form H(to) — lo/(lo + i/r). In 
the superconducting state, the frequency integral in (Q) 
can be evaluated analytically in the clean limit At 3> 1. 
After lengthy but straightforward calculations we found 
that both ai(u) and 1/t(u>) vanish at lo < 2 A, while at 
larger frequencies 
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where -E(x) is the complete elliptic integral of the second 
kind [§J. At u = 2A + 0, E = tt/2 and both ax{oj) and 
l/r(o;yjump to finite values. At high frequencies, E(x ss 
1) — » 1, <Ti(o;) vanishes as w -2 , and 1/t(lo) approaches 
the normal state result t(u>) = r. To the same order, we 
also have 11(0) = 1 — tt/(8At). We checked analytically 
that the /—sum rule (8/w 2 ;) J* +0 efcjoi (w) = 1 — 11(0) is 
indeed satisfied. 

Expanding E(x) near x — 1, we find that at high fre- 
quencies r" 1 ^) - T- 1 « (2A 2 /w 2 T)(log(2w/A) - 0.5), 
i.e. J dbj{\/T(bj) — 1/r) converges. The convergence im- 
plies that for a dirty BCS superconductor, the differential 
sum rule for l/r(w) is an exact one. 

The plots of cri(u>) and 1/t(lo) are presented in Fig. |l| 
together with the results for I a (ui) = (8/uj 2 ^ dxa\(x) 

and I t {uj) = (r/2A) dx(l/r(x) - 1/r). It follows 
from (H) that the sum rules for <7i(w) and l/r(w) should 
be exhausted at frequencies of the order of A. Fig. |l| 
shows that this is indeed the case. We next consider 
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FIG. 2: The results for 1/t(ui) (a) and <7i(w) (b) for a toy 
model in which the pairing is not accompanied by the feedback 
on the electrons. The frequency is in the units of A (the peak 
frequency in the DOS). We used A = 0.25aX Observe that 
l/r(w) is nearly homogeneously suppressed at a) = O(A), and 
the differential I t (uj) converges to zero only at u> ~ 10 3 A (not 
shown). On the other hand, 7 CT (w) converges to I a (oo) ~ 0.67 
already at u ~ 4 A. 



the behavior of ai(ui) and 1/t(u>) in the toy model in 
which we assume that the pairing does not change the 
fermionic self-energy. Having the cuprates in mind, we 
also assume, that the normal state is not a Fermi liquid, 
i.e., fermionic self-energy at low frequencies behaves as 
E(w) = (iuj) ^ 1 - 01 with a < 1, 

For the frequency dependent self-energy E(w), the 
current-current correlator II (w) is still given by Eq. (0), 
but now ft is substituted by ft + S(ft). For definitc- 
ness, we present the results for a = 1/2, which is 
the normal state quantum-critical exponent in the spin- 
fermion theory but the results are qualitatively the 
same for all a including the marginal Fermi liquid limit 
a — > |L0]]. For a = 1/2, the fermionic density of states 
N(lj) = im[£(w)/(A 2 - T, 2 (lu)) 1/2 } has a maximum at 
oj — A w A 2 /<D, that could be identified with the pseu- 
dogap. Still, N(uj) remains finite at all frequencies. The 
assumption that the self-energy does not change is in- 
ternally consistent if indeed the anomalous normal state 
behavior extends to frequencies which well exceed the 
pairing gap, i.e., when Q ^> A. 

Evaluating n(w) and substituting it into <ti(u>) and 
l/r(w) we found that in the normal state, ai jn (u>) oc 
(ljoj)^ 1 / 2 at lu <C ui and ai <n oc (w/w 3 ) 1 / 2 at uj » Co, 
while 1/T n (u)) oc (ujui) 1 ^ 2 in both limits. For A =/= 0, 
we found that for uj < A, o"i(w) oc (cja}) _1 / 2 (w/A) 5 / 2 
and l/r(w) oc (w) 1/2 (w/A) 7 / 2 . We see that both 
cri(u>) and 1/t(u) are reduced compared to their nor- 
mal state values but are still finite. At ut ~ A, both 
cr±(u>) and l/r(w) become comparable to the normal state 
results. At larger A <C u> <C Co, u\{lo) — a\^ n {uj) — 
1.992(w 2 / /47r)A(w 3 (I>)- 1 / 2 and l/r(w) = 1/t„(cj) - 

3.5lA(w/o;) 1 / 2 . Finally, at u > w, 0i(u;) — ci,n(w) oc 
w _7 / 2 logw and 1/t(lo) — l/r n (w) oc w~ 3 / 2 log w. 

We see that <j\ (to) converges to its normal state value at 
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FIG. 3: The normal state results for 1/t(w) (a) and <7i(to) 
(b) for the spin-fermion model for A = 2. The insets show 
I CT (w) and the differential I t (uj) between T = lo s j /3 and T — 
2a> 3 //3. The sum rule for 1/t(u>) is not satisfied due to weak 
convergence at high frequencies. I<j(lo) flattens at u) ~ 10a; s / 
but converges to the /—sum rule value 7 CT (oo) = 1 only at 
extremely high ~ 10 3 io a f (not shown). 



frequencies of order A, as in a dirty BCS superconductor; 
the sum rule for <j\{uj) is then exhausted at lo ~ A. This 
behavior is illustrated in Fig. |]a where we present the 
results of our numerical calculations - Ia{w) converges 
to I a {oo) = 1 - n(0) (» 0.67 for our choice of Q = 2 A) 
already at fi ~ A. On the other hand, r _1 (tj) — t~ 1 (uj) 
scales as a; -1 / 2 between u ~ A and cj ~ a) such that at 
these frequencies I r (w) = /dw(r _1 (cj) — t i 7 1 (cj)) does 
not converge. Furthermore, at these frequencies, 1/t(lo) 
is still smaller than l/r n {uS). This result holds for all 
a < 1 as one can straightforwardly verify. Only at lo > Co, 
r _1 (cj) finally becomes larger than t~ (u>), and I t (lo) 
converges. The convergence implies that the differential 
sum rule for \/t{lu) is again exactly satisfied, however 
it is exhausted only at frequencies that well exceed the 
pseudogap. We present the numerical results for 1/t{u>) 
and I t (lo) in Fig. ||b. 

We now present the results for <T\{u>) and 1/t(lo) for 
spin-fluctuation mediated d— wave pairing. We obtained 
these results by solving a set of coupled Eliashberg equa- 
tions for the spin-fermion model that describes the spin 
fluctuation exchange at low energies |D| . We will demon- 
strate that at low T, the behavior of the conductivity 
and the relaxation rate resembles that in a dirty BSC 
superconductor, while immediately below the pairing in- 
stability the system behavior is similar to that in the toy 
model for the pseudogap. 

The spin-fermion model is characterized by a single 
dimensionless coupling constant A and a single overall 
energy Q that scales with the effective spin-fermion inter- 
action ||. We will also use a characteristic energy scale 
for the spin fluctuations ui s f = lo/A\ 2 . A fit to the NMR, 
ARPES and neutron experiments yields A ~ I — 2 near 
the optimal doping ||. We refer the readers to Ref 0] 
for the discussion of the applicability of the model to the 
cuprates and the justification of the Eliashberg approach 
at strong spin-fermion coupling despite the formal ab- 
sence of the Migdal theorem. 



We begin with the normal state. In Fig. |^a we present 
our results for 1/t(lv) and I t (lj) at various T. For def- 
initeness we used A = 2. We see that I t (lo) diverges at 
high frequencies, i.e., the differential sum rule is not sat- 
isfied. We checked analytically that this is caused by the 
1/lo behavior of the integrand in I t (lo). 

For completeness, in Fig. ||b we present the results for 
ct\{uj) and I a (w) at various T. We see that io-(w) flat- 
tens up at lo > I0o; s /, but its value is still by about 
30% smaller than it should be for lo — oo. The full sum 
rule is exhausted only at unrealistically large lo ~ I0 3 w s / 
(Ref. ||) where the low-energy theory is clearly inappli- 
cable. The weak convergence of I„{lo) is related to the 
fact that over a wide frequency range o~\(lS) is inversely 
proportional to to, and I a (uS) increases as logw. The 
practical application of this result is that the plasma fre- 
quency LOpi obtained by integrating the measured ax(io) 
up to, say, leV should be about 30% less than the ac- 
tual plasma frequency. This is consistent with the fact 
that to match the measured o~i(u> — ► 0) with the theoret- 
ical <7i(w — ► 0) obtained using extracted from the 
ARPES data, one needs the plasma frequency of about 
20000cm -1 H [TTJ , while extracted from the "partial" 
sum rule |Q is around 16000cm -1 . 

We next consider what happens below the pairing in- 
stability temperature Ti ns ~ 0.2lo fll2j. Earlier Schmalian 



and the two of us have found [fL3fl that at T < Ti ns , the 
the fermionic self-energy remains large at the smallest lo 
and gradually evolves from its normal state value. This 
gradual behavior is qualitatively different from a dirty 
BCS superconductor, as in the latter the quasiparticlc 
spectral function and the density of states instantly drop 
to zero at frequencies below A due to a feedback from 
the pairing ^ Q|. We found in Jl3| that the behav- 
ior similar to that in a BSC superconductor occurs only 
below T c < Ti ns , and the difference between Ti ns and 
T c increases with increasing A. We conjectured that at 
T c < T < Ti ns , fluctuations destroy coherent supercon- 
ductivity, i.e., the system is in the pseudogap regime. 

In Fig. |i| we present the results for ai(co) and \/t(uo) 
for two different A and for three different temperatures: 
T -C T c , T c < T < T^s and T = Ti ns where <J\{lS) and 
1/t(lo) are the same as in the normal state. We clearly 
see that between T c < T < Ti ns , 1/t(u>) is nearly homo- 
geneously suppressed, while at smaller T it shifts from the 
lowest frequencies to w ~ 2A (~ 0.8a>) and develops an 
overshoot at u> > 2A. The magnitude of the overshoot 
depends on the coupling and is larger at larger A. We 
compared our results with the fermionic self-energy [jl3| , 
and found that the temperature where l/r(w) develops 
an overshoot, agrees with the onset temperature for the 
feedback effects on the fermionic propagator. 

For quantitative analysis of the sum rule we looked 
analytically into the high frequency behavior of 1/t(lo). 
We found that the temperature variation of 1/t(lv) at 
high frequencies is still strong enough to prevent the 
differential sum rule for 1/r. On the other hand, we 
see from the insert in Fig. ^b that this variation can be 
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FIG. 4: The behavior of l/r(w) and oi(u) in the spin- 
fermion model below the pseudogap temperature Ti ns for 
A = 2 (A ~ 0.3Q) and A = 1 (A ~ 0.2u>). Observe that 
the overshoot between the spectra of l/r(o>) develops only 
below T c . The inserts in the upper panels show the behav- 
ior of the differential sum rule I T (w) between T < T c and 
T c < T < Ti n3 . The differential sum rule is not exact, but 
is approximately satisfied if uj ~ 3 — 4Gj ~ 10 A. The right 
lower panels show the behavior of oi(u) vs T at small and 
large frequencies for A = 1. Observe that the changes in ai 
are confined to T c rather than to Ti na . 



safely neglected if the frequency integration is restricted 
to uj ~ 3 - Au> ~ 10A (for optimally doped YBCO, 
10A ~ 2500cm' 1 ). In Fig. |b we show the results for the 
differential sum rule between T <C T c and T c < T < Ti ns . 
We see that the spectral weight is almost exactly con- 
served if the upper limit of the frequency integral is cho- 
sen close to 3 — 4a; . The conservation of the spectral 



weight between T <C T c and T > T c , and the apparent 
loss of the spectral weight between T c and T ins are the 
main results of recent experimental analysis Q . We see 
that these results are fully reproduced in our analysis. 

We also see from Fig. [| that at small oj, a\ (u>) keeps 
increasing between Ti ns and T c . This indicates that the 
development of the pseudogap does not give rise to a 
suppression of the conductivity at the lowest frequencies. 
The latter is only reduced below T c . The sensitivity of 
o~±(u> — 0) to T c rather than to the pseudogap tempera- 
ture is consistent with the data |Q. The low frequency 
behavior of <j\ (u) = 0) below T c is not captured in our the- 
ory as it is predominantly determined by impurities |l5[ . 
Finally, we found both analytically and numerically that 
at large w, cri(w) also turns out to be sensitive to T c (se e 
the inset in Fig. ||a). This also agrees with the data p6[ . 

To conclude, in this paper we considered the differen- 
tial sum rule for the effective scattering rate 1/t(u>) (the 
difference between the area under l/r(w) for two differ- 
ent temperatures). We argued that for spin- fluctuation 
mediated pairing, this sum rule is generally not an exact 
one, but is rather well satisfied below T c and is exhausted 
at frequencies compared to the pairing gap, A. We iden- 
tified this behavior with the strong feedback from the 
pairing on the fermionic self-energy We found that in the 
pseudogap region, where feedback effects are small, the 
differential sum rule is (approximately) exhausted only 
at much larger energies, comparable to the upper limit 
of non- Fermi-liquid behavior in the normal state, while at 
u = O(A), 1/t(u>) is nearly homogeneously suppressed 
compared to the normal state. We argued that this be- 
havior as well as the behavior of <j\ (uj) are consistent with 
the experimental data on the cuprates. 
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